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Abstract 

It is shown that there exists two inner authomorpism which lead 
to different form of the sistems equations of integrable hierarchy. We 
present discrete and Backlund transformation connected with such 
systems and a general formula for multi-soliton solutions based on 
these symmetries. 
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1 Introduction 


In a series of papers by the author starting in 80’s of the previous century 
[1] a simple method for constructing integrable systems together with their 
soliton-like solutions was proposed. This method requires only two calcula- 
tional steps - solving a system of linear algebraic equations and performing 
a Gaussian decomposition of a polynomial into a product determined by its 
roots. This lead to a nonlinear symmetry of integrable system termed a 
discrete substitution or an integrable mapping. A discrete substitution is 
a nonlinear transformation (with no additional parameters) that, given an 
arbitrary solution of an integrable system, produces a new solution according 
to certain rules. 

All systems of equations invariant with respect to such a mapping are 
united into a hierarchy of integrable systems that bears the name of the cor¬ 
responding discrete substitution. This substitution is a canonical transfor¬ 
mation [2],[3] which explains successful applications of methods of canonical 
transformation theory to integrable systems [9]. 

Let us assume that an integrable system has some inner automorphism a. 
This implies that there is a solution that is also invariant under a. Usually, 
such solutions are interesting for applications. In general, a discrete trans¬ 
formation will not commute with a. However, it is possible starting from a 
non-invariant solution to produce a “good”, invariant, solution after several 
applications of the discrete transformation (implementation of this scheme is 
outlined Refs. [4] and [5]). 

On the other hand, it is well known that Backlund transformations con¬ 
tain additional parameters and therefore allow to increase the number of 
parameters in the solution after their application. Because the method of 
discrete substitution allowed to obtain n-soliton solutions directly and ex¬ 
plicitly, the author had no interest in Backlund transformations from point 
of view obtaining soliton solutions. But Backlund transformation have more 
wide sence and we would like to close this loophole in the present paper. 

The goal of the present paper is to show how to construct Backlund 
transformations for integrable systems using methods of [4],[5]. We restrict 
ourselves with the well-known example of the nonlinear Schrodinger equation 
to ease the demonstration of the main idea and methods needed for solution 
of this problem. The extension to the general case is obvious and will be 
clear after considering this single example. 
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2 Discrete transformation 


All considerations of the present paper are given for the simplest example 
of the nonlinear Schrodinger equations. However, it will become clear from 
our arguments, detailed calculations below, and references that our method 
is applicable to all integrable systems with soliton-like solutions. 

The nonlinear Schrodinger equation can be represented in two forms: the 
classical one for one complex-valued unknown function t) 

= 0 ( 1 ) 

and as a wider system of two equations for two unknown complex-valued 
functions functions u, v 

—i2ut + Uxx +‘2(uv)u = 0, i2vt + Vxx + 2,{uv)v = 0 (2) 

(The factor of 2 multiplying time derivatives can be of course eliminated by 
redehning the time variable.) The last system is obviously invariant under 
an inner automorphism cti: u ^ v,v ^ u and thus it has solutions invariant 
with respect to this change of variables. Such solutions satisfy u = v = 'ip 
and are therefore also solutions of (1). 

The system (2) is invariant with respect to the following nonlinear invert¬ 
ible transformation 

U=\ V = v{uv - {\nv)xx), ^ u = U{UV - {\nU)xx) (3) 

which is referred to as a discrete transformation, discrete substitution or an 
integrable mapping. 

Transformation (3) takes a given solution [u, v) of the system (2) to a new 
one [U, V). But this transformation is not invariant with respect to ai. Thus, 
if the initial solution has a property {u = v), the new solution {U, V) will not 
have this property. Now we make the following very important observation. 
If we apply the direct discrete transformation m times to a solution of (1) 
{u = v) with the result and apply the inverse transformation to 

the same solution m times with the result u~^, v~^, the resulting new so¬ 
lutions are related as follows: [7+”^ = u~"^, Thus, if we start 

from an obvious linear solution of the system (2) uq = 0, 2ivot + ^Oxx = 0 
and after 2m steps of direct discrete transformation obtain a solution of 
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2 m 2 mt + U2mxx = 0, U 2 m = 0 with U2m, = ^o, we are guaranteed that at the 
mth step we will obtain a (m-soliton) solution of ( 1 ) - traditional nonlin¬ 
ear Schrodinger equation. This approach together with the corresponding 
mathematical formalism is described in detail in [4],[5]. 


3 Backlund transformation 


Now we would like to hnd not the symmetry of the enlarged system (2) 
but independently the symmetry of (1). This symmetry transformation can 
contain additional numerical parameters and after each its application we 
increase the number of parameters in the solution. This is exactly Backlund’s 
original idea. However, one has to keep in mind that, applied to the general 
solution of the equation, this transformation cannot add any new independent 
parameters, but can only change the initial functions on which the general 
solution depends. 

Consider a solution of the classical nonlinear Shrodinger equation -0 = 
u = V. This equation may be written in Lax pair form 


9x9 




9t9 


/ A2 - f, 

+ -A2 + fJ 


(4) 


Now let us use the formalism of [5] and try to hnd a new solution of the 
classical nonlinear Schrodinger equation in the form 


G 


fX + A 

V C 


B 

X + D 


^9 = P 9 


(5) 


where functions A, B,C, D are determined by imposing the condition that 
matrix G has linear dependence between its columns (and rows) with con¬ 
stant coefficients ( 4 , 02 ) at two points of the complex A plane, Ai, 2 - In the 
Appendix we show that this condition is not independent but follows directly 
from (5). 

From this condition we immediately obtain DetP = (A — Ai)(A — A 2 ) 
and a linear system of equations for relating functions A, B, G, D to matrix 
elements of g, parameters Aj, and vectors d 

{Xi + A){gc)\ +B{gc)\ = {X 2 + A){gc)\ +B{gc)\ = ^ ( 6 ) 

^{ 9^)1 + (Ai -|- D){gc)l = 0, C{gc)l + (Ai -|- D){gc)l = 0 (7) 
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In particular, we have 


(Ai - A 2 ) + B{ 


{9c)l 

{ 9 c)\ { 9 c)j 




{9 c) 2 {9 c, 


^) = 0 ( 8 ) 


Now let us obtain a new L-A pair. For the matrix element {GxG ^)ii we 
have 


det 

{GxG ^)ii =- 


Ax A (A A)zA A iBv 
C 


Bx A (A A A)iu — BiX 
\ + D 


detP 


iX 

(9) 


where we used the fact that the determinants in the numerator and denom¬ 
inator have zeros at the same points. Similarly, we obtain 


det 


{GxG 


-lA 


112 


A A A 

Ax A (A A A'jiX A Biv 


B 

Bx A (A A A'jiu — BiX 


detP 


and hnally 


GxG-' 


/ A, u-2B\ 
Ua2G -a J 


i{u—2B) 

( 10 ) 

( 11 ) 


If we did not demand the condition of complex conjugation u = v* to hold, 
this would exactly be a Backlund transformation for the enlarged nonlinear 
system (2). 

But if we start from a solution of the nonlinear Schredinger equation (1) 
u = V* and would like to obtain a solution of the same equation, we have to 
require, according to the last expression, that G = —B. 

To this end, it is necessary to use known explicit expressions (8). From 
L-A representation (4) it follows that g can be considered as a unitary matrix 
gg^ = 1 under assumption that A is real. 

This means that g~^ = g^ or in a matrix form (in what follows / = /*) 


g{X)22 — g*{X)n, g{X)i2 ——g*{X)2i, 


Taking into account that the matrix elements of g are analytic functions of 
A (as solutions of a differential equation with coefficients analytic in A), we 
conclude that 


{9{^)22)* — 9{^*)lli (fi'(A)ll)* — fi'(A*)22, 
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{gW2i)* — -9{^*) i 2, {gWuT — -9{^*)2i 


Now we can evaluate C. We have 

{gcfi ^ 5'(Ai)iic} + g{\i)i2cl ^ 5^(Ai)ii + g{\i)i2ai 
( 9^)2 fi'(Ai)2ic} + g'(Ai)22C2 fi'(Ai)2i + fi'(Ai)22ai 


where ai = -f-- According to the above formulae for complex conjugation, 
we have 


(F(Ai,ai)* 


Mi 

{9c)h 


^(Ai)22 ~ g{.^*l)2,l'^l _ _ 1 

—fl'(Aj;)i2 + fi'(Aj;)iiai F(Ai,— 


For further manipulations let us rewrite (8) in terms of Fi,F 2 introduced 
above 

(Ai-A 2)+C(F(A,. ai)-F(A2. a2)) = 0, (Ai-A2)+S( = 0 

(13) 

The condition B = —C* together with the conjugation properties of functions 
F lead to 

Ai = a;, a; = (14) 

0:2 

Thus if for some solution of (1) the corresponding element g is known new 
solution of the same equation differente from the initial on the pair of complex 
parameters can be constructed by the rules of present section. (But not forget 
about the comments on this subgect in the beggining of this section). 


4 Multisoliton solutions 


Now let us apply n times the transformation of the previous section to a cer¬ 
tain solution of the nonlinear Shrodinger equation (1). Each transformation 
is dehned by two complex parameters Aj, a*. After n applications, we have 
for the corresponding matrix Gn 


G 


n 


/ A -|- A B \ „ 

V c \ + d) 


TDn 

11 

pn-1 

^21 


Pl2~^ \ 
hn 90 
-'22 / 


Here the notation P means that the coefficient at the highest degree of the 
corresponding polynomial is equal to one. Gn has zero vectors at 2n points 
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of the Xi, X* plane with the coefficients of proportionalities ctj, — For this 
reason all coefficients of polynomials of the element G„ can be obtained from 
the linear system of equations 

n—1 n—1 

A? + y +E = 0 

k=Q k={) 


n—1 

ef-\^ 


k={l 


n—1 


+ A” + i: P4A: 

k=0 


0 


where are the coefficients at in the corresponding polynomial. In con¬ 
nection with Kramers rules coefficients interesting for further considerations 
are 


■pn—l 

-'12 


dehn{xr\ A„ 1; A- F^Xr\ F,) 


pn-l 

-'^21 


def2n(A” ,..,Xi,l-,FiXi ,...FiXi,Fi) 

*fa„(A”i P-‘Ar^ ...Fr'A,, F-\ A”-‘,... A., 1) 


(15) 


(16) 


deh„(FE\r\ ...Fr‘A., F-‘i A”-',A,, 1) 

In the last formula we symbolically wrote the structure of each of 2n lines 
of the corresponding determinant matrices. In the case of a diagonal initial 
matrix qq = expi{Xx + X^t)h formula (15) was presented in [4] without any 
connection to the Backlund transformation of the present paper. 

Using absolutely the same technique as in the previous section we obtain 


2 V' “ 

Thus, after n steps of Backlund transformation each of which is dehned by pa¬ 
rameters (Aj, A*, Q(j, —we get a new solution of the nonlinear Shrodinger 
equation 

U = u- 2P^-\ V = v + 2P^G^ 

From the explicit expressions (15) and (16) it follows that all Backlund trans¬ 
formations are commutative ( functions symmetrical to per¬ 

mutation of all pairs Aj, ctj). 

If we choose the initial solution in a ’’zero” form u = v* = 0, g = 
expi2{Xx + X?t)h, we obtain an n-soliton solution in an explicit form. This 
solution of course coincides with the one obtained previously in [4],[5] using 
the method of discrete transformation and repeated in the previous section. 

The result of [4] corresponds to decomposition of the determinant of 2nth 
order into a sum of the products of the corresponding minors nth order. 
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5 Backlund transformation in its original 
form 

As it follows from the material of the previous section, Backlund transfor¬ 
mation, in contrast to the discrete one, is not local. Indeed to obtain a new 
solution U = V* from a solution u = v*, we have to resolve the L-A system 
(which is a system of linear differential equations), obtain the element g and 
use its matrix elements to construct a new solution. 

This fact is related to the original Backlund’s result, who found relations 
between the derivatives of the new and old solutions. We wonld like to show 
now, how snch relations can be obtained from the formalism of the previous 
section. 

Using the dehnition of F(Aj, ai) and taking into account equations of L-A 
pair (4), we obtain 

F: = {F{K, a),), = i{v - - u{F^f), 

F: = ^iiXsV + z|) - 2(A^ - ^)F-^ - (A.u - t^){Fr) 

Calculating F^, via 2B = u — U,2c = V — v from (13) and substitnting 
result into the system of derivatives above (with respect to a; or f arguments) 
we come to a system of two ordinary differential eqnatinos of the hrst order 
connected u, v and U, V and containing parameters Ai = A 2 in explicit form. 
Two additional parameters a!i,Q ;2 = —^ have to arised in process of inte¬ 
gration of the last system in which u = v* considered as known and U = V* 
as unknown fnnctions or visa versa. 


6 Second automorphism a 2 of nonlinear Schrodinger 
system 

Let us perform a change of variables v = e*®, u = e~^^{R — ^^) in (2), where 
R, 6 two new nnknown complex fnnction. The sence of this snbstitution will 
be clarihed something later. Instead of (2) we have 

+ (0.)' = 2R, 2Rt - - 2{e,R), = 0 (18) 
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Of course the last system is invariant with obvious exchange R ^ R*, 6 ^ 6*, 
which we call second inner automorphism a 2 of nonlinear Shrodinger system. 
The system (18) after excluding R is equivalent to single equation 

~ -^(^xxxx + ( 2 ^^ = 0 (19) 

The Author cannot say anything about the system (18) or equation (19) 
because as it seems he ever encountered them before in such a form in the 
literature. 


6.1 Method of discrete substitution 


The nature of a 2 authormorphism consists in the fact that if discrete trans¬ 
formation (3) is iterrupted on 2n -|- 1 step (but not on 2n one as it was 
in the case of Ui), it has as it conclusion (19). Indeed in this case in 
the middle of the latice arises two solutions (n„ = = '%^) and 

iun+i = = -n^)) which are connected by condition vl = Un+i = 

< = Vn+i = Vn{unVn + {\nVn)xx), which lead to (18) and (19). 

Using the technique of dicrete transformation n-soliton solution of the 
equation (19)and system (18) may be represented in the terms of the following 
notations 


2n+l 2n+l -i 

F = ^ / = ^ — 

k=l k=l 


n'(Afc - A, 


where Li = + Xx and Aj,Ci -{2n -|- 1) numerical parameters connected 

via equation F* = f. From the last condition it follows the limitation on 
parameters of the problem: 2s -|- 1 parameters = A/? are real one, the 
remaining 2(?7. — s) are in complex congugated pairs A)^ = Ab (1<A, F< 
(n — s)). In all cases CjC* = n^(Aj — A^)^. This number under all chossing of 
A parameters above is a positive number. Finally 


e 


Dn 
I In 


D. 


n+1 


R = (\n{DnDX))xx 


where Dg is determinant of the s- order of the matrix (typical for discrete 
transformation calcules): 



/ F 

F. 

\ Fxx 


Fx 

Fxa 

P 

J- XT. 


Fx 

Fx-, 


6.2 Method of the Backlund transformation 


It is necessary to repeat word by word all calculations of the section 3 up 
to (11), keeping in mind that in all formulae connected q and {u,v) now 
V = e*®, u = e-'^%R - ^0xx)- 

’’One” soliton solution, the simplest solution of (18) is the following one 
(6* = —2(AQf + Aqo;) = —2Lo, R = 0). It is obvious that correspoding group 
element g belongs to the group of lawer triangular matrices and may be 
represented in the following form: 

g{x,t : A) = 


Taking into account equations of L-A pair (4) we obtain in a consequence 

g.g-" = (ax + 2Txa)X_ + TxH = i{\H + 

gtg~^ = {at + 2Tta)X_ + TtH = i{\^H + + Ao)A_) 

Solition of the last equations are trivial with the hnally result 


1 

a = r = i{}?t + \x) 


2 A-A, 


0 


In the general case let us represent element g in ussual form of S{2,C) 
group 

' e"” + aj3e~ 


g = ^ 


ae 

o-T 


(3e-^ e-^ J 

g'g-^ = {a’ - 2r'a - P'e-^^)X+ + {R + (3'ae-^RH + pR-^^X. (21) 

Taking into acount equations of L-A pair (4) we obtain 


fA = {tx + Pxae iv = PxC 


From the last relation we conclude that pR doesn’t depend from A and 

iX — Tx 


a = 


Px^ 


-2t 
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Now let us rewrite equation defining C (8) substituting in it (20) ( to have 
no mixing we change parameters 0 : 1^2 in (12) on z/ 1 ^ 2 ) 


1 1 

7; + 1- Hit -+ «i - IT -«2) = 0 

C Ai — A2 /?i + vi P2 + ^2 


( 22 ) 


The new solution V = v + 2C (result of transformation V = v + 2C) (11) and 
condition of its inariantnes with respect to cr 2 VV* = 1 may be rewritten as 
{vv* = 1!) 


V V* 


-2 


(23) 


or using (22) we have in a consequence {iv = (3x6 ^'^) 


V (3x6 , 6^'^^ 

C i(Ai — A 2 ) (3i + z/i /?2 + ^2 


0 : 2 ) = 0 


Exept of this from (23) (n* = 1) we obtain also 


V = C(—1 + i 



1 ), 


V = C{l + i 



1 ) 


Keeping in mind all obtain above expresions and comments we rewrite 
the last equality in the form 


u I (3l 

-h 1 H-- — - 

c - \2y (3^ + pi) 






+ 1 ^ 2 ) 


- y] = 0 , 


V 1 

+^(A7^A2)^‘“(/52 + z/2)e-"^ 


[111 / . \ _ I.T 0 


Summarizing with complex congugation and not forgetting about (23) we 
come to eqution 


■[1^ 


n 




ii- A'!)! - 0 m 


i{Xi — X 2 ) {P'^ + 1 ^ 2)6 f(Ai — A 2 ) (/3^ + zz2)e "^2 


To hnd relation connected 4 up to now arbitrary numerical parameters 
Ai, 2 , 1 ^ 1,2 it is necessary to take into account second condition of conguga¬ 
tion dectated by authomrphism cr 2 , which allow to hnd connection (3 and (3* 
functions. 
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Comparing (21) and L-A pair representation (4), we obtain 


{o' — 2T'a — (3'e = iu = i-uv = i—{r — -(Int;)^;^;) 

Substituting all above obtained relations we have 

+ 2^xx = — Tx,x + + (A + iTx)9x 

In a L-A pair formalism A have to be considered as a real parameter. Thus 
imiginary part of last equation has the form 


Oxx = -(r* - t)xx + {t* + t)x9x + {Tx 


? - 


The last equation has obvious first integral 


Ox - -{t* - t)x = 


Reminding that ie^ = jSxC and substituting into the last quation we trans¬ 
form it to the form 


1 


-i}^Pxe 


-(r*+r)y^ _ 1 f/9„p-C*+'r)V _ 


^ I 

Second integration leads to hnal result 
1 


= ciPx 


-{(3xe-^-*+^)) = = Ci/3 + Cq 


(25) 


we pay attention of the reader that in (25)parameters Ci, Cq not depend from 
(x, t) variables, but only from real parameter A. Inverting (25) we obtain 
(3 in form in which it became clear its structure with respect to complex 
conjuguation 


1 


+ z/= + z/-— 

Cl Cl 

Putting in the last expression indexes /c = 1,2 for A we come to elements 
introducing in (24) 

I3k + i^k = + h. 

Cl 
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where Pk = i^k - 


(^l + ^l = 






After substitution the last expressions into (24) and simplest manipulations 
we come to two possibilities. First one Ai = AT, A 2 = AS, hihf = i/iw , ^ 2^2 — 

and the second one Ai = A 2 , hii >2 = Of course this result in the 

case of multisoliton solutions coinsedes with obtained above by the method 
of discrete substitution. 


7 Outlook 

First, we would like to explain why, in spite of our comments in the Intro¬ 
duction, Backlund transformations are interesting objects for investigation. 
Backlund transformations, as well as discrete ones, are canonical transforma¬ 
tions. This means that under these transformations densities of conserved 
quantities change by complete derivatives. In the case we considered above 
this applies in particular to the energy density. Therefore, if the energy 
density of the initial solution is e, after the application of the Backlund 
transformation it will be = e -f A^,. The term A contains all parameters 
Ai = A 2 , which dehne the Backlund transformation. For obtaining 

an explicit expression for A an additional (not a very cumbersome) calcula¬ 
tion, using the technique of [4],[5], is necessary. If e itself has been obtained 
from some other solution and so on, we will have a hnal expression 

Each A* contains only parameters of previous transformations. Because the 
total energy for an n-soliton solution is equal to n, it is very natural to 
assume that all terms of the last sum have the same behavior at inhnity and 
contribute 1 to the total energy after integration. 

It is completely obvious that ai and a 2 exist for all integrable systems 
possessing soliton solutions and a discrete transformation. Indeed, as it was 
explained above, this is related only to the parity of the step, odd (2n) or even 
(2n-|-l), on which the discrete transformation applied to the initial solution 
is interrupted. 
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Without any doubts all results of this paper also apply to supersymmetric 
integrable hierarchies [ 6 ] and multicomponent matrix-type hierarchies [7],[ 8 ]. 

What are the systems connected to (72 is completely unknown to the 
author at this moment. 
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In this Appendix we would like to show that 
Let us rewrite (4) and (5) in equivalent form 


Bx 

Cx Dx 


+ 


A -h A 
C 


B 

\ + D 


A, 

V 


u 

-A 


= i 


A, U 
V -A 


A +A 
C 


B 

\ + D 


At Bt 
Ct Dt 


+ 


A +A B 
C X + D 


A^-f, An-if 
An + i^ -A2 + ^ 




XU-t' 

XV + i% -A2 + ^ 


A A 
C 


2 ^ 


2 ''' 2 

From which we have the following system of equations 


B 

X + D 


u—U = 2B, v—V = —2(7, iAx = Bv—CU, iDx = Cu—BV, iBx = Au—DU, iCx 


and the same kind of the equations with respect to t differetiation. As result 
we have A -|- D,AD — BC doesn’t depend from x,t arguments and thus 
representable in the form A -|- D = Ai -|- A 2 , AD — BC = A 1 A 2 . Besides this 
functions (7, D satisfy the following system of equations 

-iCt - ]^{Cxx + ^tkCx + Ak^C) - C{BC) + = 0 

iBt - ^{Bxx - UkBx + Ak^B) - B{BC) + = 0 


Dv—A 
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where A: = Ai + A 2 , e = 

The last system of course is equivalent to nonlinear Schrodinger equation. 
More other solution of (2) is connected with the solution of the above system 
via relations 


v = -C- 


iCx + kC 
2^/e-BC' 


u = B + 


iBx — kC 
2Ve - BC 
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